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Abstract 

In the paper [V. Adler, IMRN 1 (1998) 1-4] a lattice version of the Krichever-Novikov 
equation was constructed. We present in this note its Lax pair and discuss its elliptic 
form. 
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1 Introduction 



Integrable nonlinear evolution equations that can be solved in terms of linear problems 
with spectral parameter on an elliptic curve have been known since a number of years, the 
most well-known example being the Landau-Lifschitz (LL) equation, cf. [||], i.e. 

St = S X S^^ + S X JS . (1.1) 

Here S is the normalised spin- vector S G R^, and J = diag( Ji, J2, J3) is the diagonal 
matrix of anisotropy parameters. 

Another well-known example is the Krichever-Novikov (KN) equation which reads 

6 = y... + ^^^^^ - ^pme. , (1-2) 

with p being the Weierstrass p-function. Eq. ( |1.2| ) arose in Q from the general problem 
of describing finite- gap solutions of the Kadomtsev-Petviashvili (KP) equation associated 
with commuting differential operators in the KP hierarchy. Eq. ( |1.2D applies to the partic- 
ular case that these operators describe a subalgebra generated by second order differential 
operators, i.e. the common eigenspaces are of dimension 2 (i.e. the "rank 2" situation) 
and obey an algebraic relation which determines an elliptic curve, cf. also [^. 

In its equivalent rational form the equation reads 

1 3 r(u) - ui^ ,^ 
ut = -,u^xx^^— 1.3 

for u = p(^) and where r(u) = Au"^ — g2U — is the polynomial of the standard Weierstrass 
curve 

r : C/2 = r{u) = 47.3 - g^u - 53 (1-4) 

(in principle by performing homographic transformations of the type 

au + P 
u 

'JU + 

r{u) can be replaced by a general quartic polynomial ). Degenerate cases (when the curve 



(1.4) reduces to a rational curve) can be mapped to the Schwarzian KdV (SKdV) equation 



= {z,x] , {z,x}^^--^ , (1.5) 

^ ^X ^ ^x 
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the brackets denoting the Schwarzian derivative. The latter equation has a beautiful 
discrete counterpart on the 2-dimensional lattice, namely 

i^n,m Zn-\-i {Zfi rn+l ^n+l,m+l) Q /-, „\ 

7 Z v Z T ^ ~ ^ 

which (together with its Lax pair) was first given in Q, cf. also The r.h.s. of eq. 



(1.6) taking the form of the canonical cross-ratio of four points in the complex plane, this 
discrete equation turned out to play the defining role in a theory of discrete conformal 
maps, cf. [0]. It also played an earlier role in the theory of Pade Approximants, cf. e.g. 



As far as integrable discrete systems associated with elliptic curves is concerned, a 
number of results exists on nonlinear chains, i.e. differential-difference equations (DAE's), 
mostly due to the group from Ufa, cf. e.g. One particular case of such a DAE is 

the following example due to Yamilov, cf. 

^ = {iul - g,u^ - 53) ( + ) , (1.7) 

at \Un+l-Un Un-Un-lJ 

which constitutes a differential-difference analogue of the KN equation. The elliptic form 
of this equation is given by: 

dtin = Q{in + in+l) + C(^n " Cn+l) " Q^n + ^n-l) " C(Cn " ^n-l) (1-8) 

with Un = p{£,n) and where C denotes the Weierstrass ("-function. In this context an ex- 
tension of the theory developed in Q to commuting difference equations that was recently 



proposed in |10| should be mentioned as well. 



On the level of integrable partial difference equations (PAE's), i.e. equations with 
both spatial as well as temporal independent variable discrete, the situation is much less 
developed. To my knowledge the first example of a fully discrete system (discrete in 



space as well as time) was given in |11|, where a lattice version of the LL equations was 
constructed from an Ansatz for a Lax pair. This resulting lattice system can be retrieved 
from the basic equations: 

TqS-S xJT = TqS -jfxS (1.9a) 
SqT-TxJS = SqS -JSxT (1.9b) 
SoTo-foSo = a(s ■ J-^T -T ■ J-^s) (1.9c) 



with a = J1J2J3 and where Sq, Tq are scalar fields and the three- vectors S, T in are 
normalised to have unit length. How this system is related to other discretisations of the 
LL equations, cf. e.g. ||l^, is yet unknown. 

In the present letter we investigate a more recent and simpler example of a PAE 
associated with elliptic curves that was found by V. Adler in on the basis of Backlund 
transformations (BT's) for the KN equation. The main point of this letter is to present 
the Lax pair for this lattice system and to exhibit its connection to other lattice systems, 
notably its elliptic form. 



2 The lattice KN System 



A lattice version of the KN equation ( |1.3| ) was found recently by V. Adler in |13| as the 
permutability condition for the BT's of the KN equation. For the purpose of fixing the 
notations, we we briefiy recall here its construction. 



The form of the auto-BT mapping solutions of eq. ( |1.2D to solutions of itself was 
proposed to be of the following form: 

UxUx = ^H{u,u,a) (2.1) 

with 

52 ^ 2 



H{x,y,z) = i^xy + xz + yz + —j - {x + y + z){4xyz - gs) (2.2) 
and where (a. A) = (p(a), p'{a)) is a point on the elliptic curve F. 



Remark: As indicated in |13] this result can be generalised to include non-auto- 
Backlund transformations between solutions of different versions of the KN equations 
associated with different elliptic curves. In fact, if h{u,v) is any biquadratic polynomial, 
then the equation h = defines a birational correspondence between the curves 10"^ = r{u) 
and = R{v) obtained from: 

r(n) = hi — 2hhvv and R{v) = ~ '^hhuu ■ 
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The permutability condition for the BT, i.e. the condition that the BT given by 

UxUx = -^H{u,u,b) (2.3) 

(where again (6, B) = {p{(3), p'{(3)) is a point on the curve T) commutes with the original 
BT leads to the algebraic equation of the form: 

kQUuuu — ki (^uuu + uuu + uuu + uu'uj + /c2 (uu + wij 

—ks (uu + uuj — ki iuu + uu \ + k^iu + u + u + u] + fee = . (2-4) 



In fact, eliminating the derivatives from ( |2.lD and (|2.3| ) one obtains the expression: 

B^H{u, u, a)H(u, u, a) - A^H{u, u, b)H{u, u,b) = 

the l.h.s. of which factorises into two branches, one being the l.h.s. of ( p.4| ) and the other 
being a similar expression with B —B. Thus one obtains the parametrisation of the 
coefficients: 

ko = A + B , ki = aB + bA , k2 = a^B + b^A 

2 

^5 = TT^o + , fee = ^ko + gski 

2 4 lb 

where both {a, A) and {b,B) are points of F, i.e. we have the relations A"^ = r{a), 

B^ = r{b) . 

In eq. ( |2.4| ) we have adopted our preferred short-hand notation for lattice systems 
with u = Un,m, u = Un+i,m, u = Un,m+i and u = Uri+i,m+i denoting the values of the 
dependent variable u around an elementary plaquette, cf. Figure 1. The a, (3 in Figure 1 
are lattice parameters representing the grid size. Alternatively, we can think of the lattice 



parameters in the case of ( |2.4|) to take values as points on the elliptic curve (1^), i.e. they 
are parametrised by a and j3 as 

(a, A) = (p(a), p'(a)) , (6, B) = (p(/3), p'(/3)) . 
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Figure 1: Configuration of lattice points in the lattice equation of KdV type. 

We mention tliat Adler's lattice system (^^) can be written in the following remarkable 
form: 

A [(n -b){u-h)-{a- b){c - b)] \{u - - 6) - (a - 6)(c - b) 
+B [{u — a){u — a) — {b — a)(c — a)] {u — a){u — a) — {b — a)(c — a) 
= ABC{a - b) (2.6) 

in which c = p{fj — a), C = p' {j3 — a). To obtain ( |2.6D we note the following relations 
among the coefficients ki = ki{a,f3) in (|2.5|) 



as well as 



ks + k2 = B{a - b){a - c) , k4^ + k2 = A{b - a){b - c) 



k5 = a{k2 + k-i) + b{k2 + k^) - Ab^ - Ba^ , 
fee = A{b^ -{a-b){c-b)f + B{a^ -{b-a){c-a)f 
+AB[A{b-c) + B{a-c)] , 

eliminating the moduli §2, gs- Use is made of the addition formulae for the Weierstrass 
p-function in the form: 

A{c-b) + B{c-a) = C{a-b) . 



We explain in the next section in what sense we consider an equation of the type (2.4) 
to be integrable, and present a derivation of its Lax pair. 



3 Integrability of the Adler system 



We will here recollect a few general points regarding lattice systems of the type presented in 
the previous section. The points made here are not new and were presented and illustrated 
for different lattice systems at numerous occasions, cf. e.g. |14, 



Eq. (|2.4| ) is an example in "general position" of a large class of partial difference 
equations (PAE's) of the following canonical form: 

f{u,u,u,u-a,(3) = 0, (3-1) 

where we adopt the notation of the previous section to indicate the vertices around an 
elementary plaquette on a rectangular lattice: 

U ■ — Un,m ) U — tin+l^m 
U '■= Un,m+1 , U = Un+l,m+l 

cf. Figure 1. In eq. (|3.lD a, (3 denote (as in the previous section) lattice parameters, i.e. 



parameters that are singled out from other (fixed) parameters by being associated with 
the lattice shifts 



/9 ^ 
U — > U 



In fact, we will assume in what follows that the lattice equations (3T) are covariant with 
respect to the simultaneous interchangement of both lattice variables as well as lattice 
parameters. As is well-known integrable lattice equations of the type above can be reinter- 
preted as the permutability condition of Backlund transformations, identifying the lattice 
parameters as the coresponding Backlund parameters. Furthermore, the dependence on 
these parameters can be exploited to obtain continuum limits of the lattice equations to 
semi-continuous and fully continuous equations that are compatible with the original fully 
discrete equation, cf. [l8| . 

Examples of ( |3.lD have been studied since many years, cf. for a review. Two 
important questions that arise in considering integrable cases of such equations are the 
following: 



1. When does a PAE of the form ( t3.l| ) make sense from the point of view of initial 
value problems (IVP's)? This is the question on the well-posedness of the IVP on 
the lattice. 

2. What is the characteristic property of such equations associated with its integrabil- 
ity? We have in mind here the analogous property to the existence of conservation 
laws and higher symmetries. 

Both questions have been answered (for canonical examples) in recent years. The first 
question was answered in a series of papers a decade ago, cf. e.g. |2^. Clearly, if the 
equation / = can be solved uniquely at each vertex of the plaquette (i.e. if it is linear 
in each of the dependent variables associated with its vertices), we can define initial value 
problems on configurations like: 




Figure 2: Configuration of initial data on the lattice for lattice equations of KdV type. 

The second question, regarding the characteristic of integrability of the lattice equa- 
tions, was discussed in (as well as in various public lectures presented by the author 
over the past few years). There, the main characteristic of the integrability is expressed 
by the fact that in the integrable cases the equation ( [j.lj ) represents actually a parameter- 
family of consistent equations on a multi- dimensional lattice. It is at this point that the 
lattice parameters a, (3 will play a crucial role: instead of fixing them once an forever and 
solving the equation for specific values of the parameters, we leave them free and associate 
new lattice directions with each new value of the corresponding lattice parameter. The 
fact that this can actually be done in a consistent way, imposing a multitude of equations 
of the same form (but with different parameters and in different latice directions) on one 
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and the same dependent variable u is the true characteristic of the integrabihty of the 
lattice equation: it is a purely combinatorial property of the lattice system that forms the 
precise analogue of the existence of higher symmetries of continuous evolution equations 
of KdV type. 



Proposition: The lattice equation (2.4), with the parametrisation of coefficients given 



by (|2.5D , represents a compatible parameter-family of partial difference equations, i.e. the 
equation can be embedded in a consistent way into a multidimensional lattice, on each 
twodimensional sublattice of shifts of which an equation of the form ( p.4| ) holds for one 
and the same independent variable u. 



Let us explain how the consistency of this embedding manifests itself, cf. |14]. As 
stated above, with each value of the lattice parameters a, /3 we can associate a discrete 
variable corresponding to a direction in a multidimensional lattice, s.t. the solution u can 
be considered as a function on this multidimensional lattice, i.e. a single-valued object 

u = Un,m,h,... = u{n, m, /i, . . . ; a, 7, . . . ) 

with the notation for the lattice shifts 

U . — ^n+l,m,fc,... ) U . — 'Un,m+l,h,... j ■ — '^n,m,h+l,... 

It suffices to investigate how the embedding the lattice equation works in a three- 
dimensional lattice. Thus we impose the system of three equations: 

f{u,u,u,u;a,(3)=0 (3.2a) 
f{u,u,u,u;a,^)=0 (3.2b) 
f{u,u,u,u;P,j) =0 , (3.2c) 

and investigate the iteration scheme around an elementary cube for given initial values 
at u, u, u and u (see Figure 3). Iterating along the faces of the cube by using the three 
equations (|3^ ), (|3?2b| ) and ( |33^ we arrive at a possible point of conflict for the value 



of u which, in principle, can be calculated in three different ways. If, generically, these 
three different ways give the same answer, i.e. if the evaluation of the iteration scheme at 
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that corner of the cube is single- valued, the system is consistent 
consider the original lattice equation to be integrable. 



It is in this case that we 
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/ ' u 
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Figure 3: Consistency of the lattice system ^-4 ) on the three-dimensional lattice. 



For all previously considered lattice equations, e.g. the case ( jl.q ), the consistency in the 
sense explained above can be easily verified by hand, and calculating the variable at the 
outer point of the elementary cube in terms of the initial values leads to an expression 
that is invariant under the interchange of lattice directions, cf. [21| for details, (in the 
paper we gave the lattice modified KdV equation as an example, where the scheme 



works precisely in the same way). We mention that in a recent submission, |22], the 
consistency around the cube was reiterated in a slightly different context. The relevant 
formulae pertaining to the case of the lattice Schwarzian KdV of , eq. (|1.6| ) , were given 
in the thesis [^], and have also been recovered in that paper. 

Coming back to the main issue of the present paper, the integrability of the Adler 



system (2.4), in the sense of the consistency on the multidimensional lattice, is straight- 
forward but technically much more involved than the previously considered cases of KdV 
type. The difficulty stems from the algebraic relations between the various lattice pa- 
rameters on the elliptic curve, and consequently the use of algebraic manipulation tools 



becomes inevitable in this case. The verification that the Adler system (2.4) obeys the 
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consistency around the cube has been verified exphcitely with the help of J. Hietarinta 
using REDUCE. Thus, in the sense of this notion of the existence of commuting discrete 
flows in the multi-dimensional lattice, the Adler system is integrable. 



We will now derive the Lax pair for the Adler system from the above considerations. 
This is achieved by noting that, as an immediate consequence of the statement that the 
integrable lattice equations are consistently embedded in a multidimensional lattice, one 
can state that in a sense (viewed as a parameter-family of equations rather than a single 
equation on a 2D lattice) the lattice equation forms its own Lax pair. This point of view was 
presented in a number of recent seminars on the subject, cf. and was implemented 



explicitly for various cases (including e.g. the case of the lattice Boussinesq equation), cf. 
also the thesis where the derivation for the case of the lattice SKdV equation ( |1.6| ) 
was presented. Let us outline the steps needed in order to implement this idea; they are 
as follows: 

1. Fix a "virtual" direction on the 3D lattice, e.g. the direction associated with a shift 
denoted by u — > u and a lattice parameter k, with corresponding equations given 
by (l33b| ) and ( ^ ). 



2. Regard now the shifted object in the virtual direction n to be a new quantity, while 
u and u represent the "physical" shifts on the original variable, and solve for u and 
u (which can be done due to the multilinearity of the lattice equation). This yields 
expressions which are fractional linear in u. 

3. Resolve the resulting "discrete Riccati" equations by linearisation, i.e. set u = f /g 
and separate into linear eqs. for / and g, leading to a 2x2 matrix system; this will 
produce the Lax matrices up to a common factor D which needs to be specified by 
a determinantal condition. 

In this way we obtain in a straightforward way the Lax pair for the given lattice equation 
in the form 

m = Lip 

(3.3) 

if = M(f 
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with (/9 = (/, g)'^ , and where the auxihary lattice parameter k plays now the role of the 
spectral parameter. Of course, the above scheme only applies to the case where the lattice 
equation can be consistently embedded into a multidimensional lattice. 



Thus, having assessed the integrability in the above sense for the Adler system 
it is an easy exercise to implement the steps 1-3 for this lattice system. This leads to a 
Lax representation of the form (|3.3D where one part of the Lax pair is given by 



if = Li^{u, u; a)ip (3.4a) 

with Lax matrix: 



1 / k-iuu — kou + kdU — , k-iuu — k^(u -\- u) — kp. 
' kouu — ki{u + u) — k^ , —kiuu + k2U — k^u + k^ 



The coefficients in ( |3.4bD can be inferred immediately from the form of the lattice equations 
itself, setting ki = ki{a, k) , i.e. replacing (3 by the spectral parameter k in the formulae 
( |2.5D for the ki. The only subtlety is the choice of the prefactor D = D{u,u,a, k) which 
follows from the determinant of the matrix in ( 3.4b| ), which is calculated to be: 



{a - kfKK'H{u,u,a) with k = p{k) , K = p'{k,) , K' = p'{K,-a). 

Choosing D to be equal to the square root of this expression to ensure that det(LK) = 1- 
Taking the other part of the Lax pair in the same form apart from the obvious replace- 
ments: ~ ^ ^, a — > /?, we have the compatibility relations 

LK(n, n; a)Lf^{u, u; (3) = L^{u, u; P)L^{u, u; a) 

which are the Lax equations for Adler system (p.4|). 

One further exercise is to obtain the Lax pair for the BT (^]^) itself. Obviously, we 
can think of this equation as a differential-difference equation representing a continuous 
commuting flow to the fully discrete lattice system ( ^.4| ) . The corresponding Lax matrix is 
obtained in precisely the same way as the discrete Lax matrices: solve the Riccati equation 
associated with the virtual shift ~, i.e. eq. (|2.1| ) with a replaced by k, A by K, and u by 
u, and then making the same identifications as before to obtain a 2x2 differential matrix 
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system. There is now an additive freedom built into the system, which is fixed by imposing 
that the Lax matrix be traceless. Thus, we obtain an additional linear equation given by 

= Ui^ip (3.5a) 

with 

1 / \9:i-{u + k){uk-\g2) gz{u + k) + {uk + \g2f 
Uk = — (3-5b) 
\ -{u - kf -i<73 + {u + k){uk - \g2) 



which supplements the discrete Lax pair (3^). The linear equation (3^) is the spatial 



part of the Lax pair for the continuous KN equation (|1.3|), which can be recovered from 
the original Lax pair given in the paper 



Thus, we have obtained the Lax pair the lattice KN equation ( |2.4D of Adler directly 
from the equation itself. After having obtained these results, which were presented during 
the workshop on Discrete Systems and Integrability at the Isaac Newton Institute (Septem- 



ber 2001), cf. |16], V. Adler communicated to the author that similar (yet unpublished) 



results had been obtained by Yu. Suris and himself. 



4 Discussion 

In this paper we presented a novel Lax pair for the lattice Krichener-Novikov equation 
of Adler. The derivation is based on the simple observation that in a precise sense the 
lattice equation is its own Lax pair. This is based on the nontrivial fact that the lattice 
equation under consideration represent a consistent parameter-family of equations which 
are embedded in a lattice of arbitrary dimensionality. The "consistency around a cube" , 



cf. Figure 3, which demonstrates this property was first presented in a paper [14| and 
the derivations of Lax pair in cases such as the lattice BSQ and SKdV equations was 
implemented in the thesis [^. In the recent submission some of these results were 



given a geometrical context, equation that were already presented in |21|. 

Although the above observation is very simple, it is also quite deep in that it consti- 
tutes really the precise anlogue of the existence of hierarchies and of infinite sequences of 
conservation laws that are the hallmark of soliton systems. 
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Having obtained the Lax pair for the Adler system one can now embark on a more 
systematic study of this new elhptic system, including its finite-dimensional and similarity 
reductions leading to finite-dimensional mappings associated with spectral problems on the 
torus and (possibly) to Painleve-type of equations possessing isomonodromic deformation 
problems on the torus. The first part of such a study is already underway, cf. [ p^ . 

Let us finish by making some remarks on the elliptic version of some of the equations 



mentioned. It is easily noted that the BT {2A) from which the Adler system originated 
can be cast into elliptic form by going over to dependent variables u = p{(,) as in the case 
of the continuous KN equation (i.e. the correspondence between the forms ( p..2| ) and (|1.3|). 
We, thus, obtain for ^ the equation: 



(4.1) 



This elliptic form can be found from the original equation for the BT by employing a 
number of relations, notably the identity: 



H{u, V, a) = {u — vY 



I fU -V 
4 V u — t> 



{u + v + a) 



1 fu + v 

4 



u — V 



{u + v + a) 



(4.2) 



in which C/^ = r(ti), V"^ = r{v)^ as well as the addition formula for the Weierstrass 
p-function, i.e. 

2 



p(0 + piv) + + v) 



1 f p'iO-p'iv) 
piO - piv) 



and the relations 



P(0 - P{r]) 



cr{r] + C)cr{r] - C) 



in which a is the Weierstrass u-function. 



p'iO 



am 



Eliminating the derivatives from the BT (|4.1| ) together with a second BT with 
replaced with ^ and a by /3 we obtain 

ct(^ + ^+ a)a{i + I- a)a{i - ^+ a)a{i - ^ - a) 



o"(^+ C + a)a(C + C- a)a(C- ^+ a)a{^ - ^- a) _ a^{2a) 



^(e + C + PM^ + ^ - PM^ - e + PM^ -i-P) 



a\2i3) 



(4.3) 
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We remark that eq. (|4.3|) cannot yet be regarded as the elhptic form of the Adler system 



(2.4), as the latter is obtained from a factorised form of the permutabihty condition of 
BT's. In fact, an elhptic version can be obtained in a straightforward way from the 
formula (|2.6| ), but we haven't managed yet to cast it into an elegant form, so we omit it 
here. Alternatively, a so-called "three-leg" form of the Adler system was announced in 
the recent submission [^2|, but it is yet unclear how this is related to the formula ([4.3|). 
The purpose of these ways of rewriting the Adler system is to obtain a better insight 
into the underlying structure of the equation, including the origin of the Lax pair. A 
more systematic understanding of the relationship between the various members in the 
KN family of equations seems to be necessary and we intend to address these problems 
more at length in a future publication, [25|. 
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